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Abstract 



We derive an analytic, albeit approximate, expression for the ground state en- 
ergy of N Coulomb interacting anyons with fractional statistics v, < \v\ < 1, 
confined in a two-dimensional well (with characteristic frequency ujq) and sub- 
jected to an external magnetic field (with cyclotron frequency u c ). We apply a 
variational principle combined with a regularization procedure which consists 
of fitting a cut-off parameter to existing exact analytical results in the non- 
interacting case, and to numerical calculations for electrons in quantum dots 
in the interacting case. The resulting expression depends upon parameters 
of the system N, uq, ro, as and u c , where ro represents a characteristic 
unit length and a# the Bohr radius. Validity of the result is critically as- 
sessed by comparison with exact, approximate, and numerical results from 
the literature. 
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I. INTRODUCTION 



The exchange statistics of particles whose orbital motion is restricted to two space dimen- 
sions (2D) differs substantially from the 3D case. The topology of their (multiply connected) 
configuration space allows for fractional statistics [1] , characterized by a continuous param- 
eter v which labels the possible one-dimensional representations of the braid group. For 
particles in 2D, v may attain values between (for bosons) and 1 (for fermions), thus 2D 
particles are called anyons [2]. The concept of anyons has been used to describe quasi-particle 
excitations in the fractional quantum Hall regime [3-5] and in high-T c superconductors [6]. 

A system of particular current interest is that of 2D electrons in a parabolic confinement 
potential, the so-called quantum dot or artificial atom [7]. These systems, realized in semi- 
conductor nanostructures, are objects of fundamental studies of ground state properties of 
interacting iV-particle systems and have also a potential for applications in quantum infor- 
mation and computation [8] . Exact closed-form solutions of the problem is reduced to a few 
simple cases due to its intrinsic mathematical complexity [9,10]. Typically, ground state cal- 
culations make use of numerical simulations for individual choices of the parameters of the 
system. It would be desirable to have an accurate, albeit approximate, analytical expression 
describing the ground state energy of such a system as a function of the parameters of the 
system (including an external magnetic field applied perpendicular to the plane of the dot). 
In order to derive such a formula we make use of the anyon concept including the effect of 
the Coulomb interaction. 

Anyons in a parabolic confining potential (with and without an external magnetic field) 
has been the subject of several investigations in the past. For the N = 2 case an exact 
solution to the spectral problem exists [1-3,5]. For the noninteracting case, its generalization 
to iV = 3 (without magnetic field) is considered in [11]. Numerical calculations have been 
performed for the lower part of the energy spectrum for N = 3 [12] and N = 4 [13]. 
Making use of the separation of the center-of-mass and relative motion one finds in the 
two-anyon case that the ground state energy for the relative motion is a linear function of v 
for < v < 1. For N = 3 and 4 [12,13,5] the ground state energy of the relative motion is 
a linear function of v only near the bosonic limit v ~ 0. The ground state in the fermionic 
end (u = 1) is continuously connected to an excited state of the bosonic spectrum (u = 0) 
and consequently - when adding the lowest energy of the center-of-mass motion (which does 
not depend on v) - one finds for N = 3,4 (and likely for all higher N [14]) a nonlinear v 
dependence of the ground state energy. In the case of interacting anyons in no external 
magnetic field the energy spectrum for two anyons at some fixed values of the Coulomb 
interaction parameter was found analytically in [15] and approximately, for two and three 
anyons, in [16]. 

In the presence of an external magnetic field one can study the interplay between the 
statistical and the physical magnetic flux. This has been done analytically for the N=2 
noninteracting anyon case with confinement [17], while the cases with N > 2 have been 
investigated preferentially without confinement [18-20]. The ground state for the N = 3 
case including confinement and magnetic field was calculated in Ref. [21]. The case with 
the applied external magnetic field and Coulomb interaction for two anyons in a harmonic 
potential was considered in Ref. [22] and for two and three anyons in Ref. [23]. Ref. [5] 
provides a review about all these studies. 
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In our treatment of the iV-anyon problem we make use of the bosonic representation 
of anyons that works with a gauge vector potential to account for the fractional exchange 
statistics but allows to use a product ansatz for the iV-body wave function. We apply a 
variational principle by constructing this wave function from single-particle gaussians of 
variable shape. It is well-known from perturbative ground state calculations for anyons in 
an oscillator potential that the expression for the ground state energy has a logarithmic 
divergence connected with a cut-off parameter for the interparticle distance [24-28,14]. We 
face the same problem in our variational treatment. Making use of the physical argument 
(see Ref. [4]) that for v ^ this distance has to have some finite value, we regularize the 
formula obtained for the ground state energy by an appropriate procedure that takes into 
account some existing exact analytical results, in the case without Coulomb interaction, 
and numerical results for electrons in quantum dots, in the case with Coulomb interaction. 
Our formula, which is an approximate closed-form expression depending upon u, N, u 
(confinement parameter), r /a B (Coulomb interaction parameter), where r = (h/ (Mujq)) 1 ! 2 
and as = h 2 /(Me 2 ) (in the presence of a magnetic field also of the parameter oj c /ujq, and 
we need to replace v by \u\), will be compared to exact, approximate and numerical results 
for quantum dots reported in the literature. 

The paper is organized as follows: In Section II we describe the system and motivate the 
ansatz for the variational treatment, in Section III we present the calculations without, and 
in Section IV with a homogeneous magnetic field for the case without Coulomb interaction. 
Calculations including Coulomb interactions are presented in Sections V and VI. Finally, 
Section VII summarizes the main conclusions. 

We would like to note the existence of two seemingly unrelated notions of anyonic statis- 
tics in the literature. One originally introduced in first quantization in the coordinate 
representation, and another derived within the framework of quantum field theory. In both 
cases the original motivation to introduce such particles was basically as an inherent possi- 
bility in the kinematics of (2+l)-dimensional quantum mechanics and clearly the concepts, 
if correctly implemented, should be equivalent whether one uses first quantization in the 
coordinate representation or second quantization [29,30]. Within the framework of quantum 
field theory fermions can be kinematically transformed into hard-core bosons (through sta- 
tistical transmutation) but not into canonical ones, thus preserving the exclusion statistics 
properties of the particles. More generically, the Hamiltonian spectra of particles sharing 
the same exclusion statistics can be connected through a continuous mapping. The anyon 
notion used in the present manuscript is consistent with the one developed in the framework 
of the bosonic representation in first quantization. Had we used the fermionic representation 
we would have ended up in an excited bosonic state. 

II. INTERACTING ANYONS IN A 2D PARABOLIC WELL IN THE PRESENCE 
OF AN EXTERNAL MAGNETIC FIELD: GENERAL SETUP 

The Hamiltonian of iV spinless anyons of mass M and charge e confined to a 2D parabolic 
well, interacting through Coulomb repulsions, and in the presence of an external homoge- 
neous magnetic field, H = He z {uj c = \eH\/(Mc)), is given by 
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1 N r 2 i 1 N p 2 

H = ^H {Pk- (Mn) + eA ext (n)/c)) + M^l\f k \* +-£-_. (1) 

k,j=ik I I 



fe=i 



Here and j4 represent the position and momentum operators of the fcth anyon in two 
space dimensions, 



N 



^)=^El^ (2) 

is the anyon gauge vector potential [11,31], r^j = rk — rj, and e z is the unit vector normal 
to the 2D plane. The factor v determines the fractional statistics (or spin) of the anyon: it 
varies between v — (bosons) and v — 1 (fermions).The external magnetic field enters by 
minimally coupling the vector potential A ext (fk) = H x fk/2. 

In order to find an analytic expression for the ground state energy as a function of is, N, 
u>o, r /a B and uj c /ujq (in the presence of a magnetic field v is replaced by \v\ (see Section IV 
below)) we employ a variational scheme by minimizing the expression for the total energy 

E = J** t (R)H*t(R) dR 
fy* T (R)V T (R) dR ' 

with a trial wave function ty T (R) depending on the configuration R = {ri....r N } of the iV 

— * 

anyons. To motivate the choice of ^t(R) we invoke the mean- field approximation to the 
gauge vector field 

X(r) = ^B u xr (4) 
introduced by Fetter, Hanna, and Laughlin [32]. This single-particle vector potential can 

— * 

be understood as that of a homogeneous "magnetic" field B u = 2nphh , e z connected with 
the carrier density p and the anyonic factor v (note: B u vanishes in the bosonic limit). By 
analogy to a physical magnetic field one can introduce a "magnetic" length l v = (h/ 'By) 1 / 2 . 
The other characteristic length of the system is the mean distance between particles r = 
1 / y/iTp. Taking into account only this mean gauge vector field (and not the external parabolic 
confining potential) one obtains a Landau spectrum [33] and it is reasonable, in the bosonic 
representation of anyons when the many-body wave function takes the product form 

N 

* T (R) = II Mr k ) , (5) 
k=i 

to adopt the single-particle trial functions iprifk) in the form 

Mn) = Cexp (-(a> + v) ^ 1 ^ 1 ^ (6) 

typical for the lowest Landau level. Here C is a normalization constant and a' a variational 
parameter. To include the external confining potential we identify r with the characteristic 
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length {h/ Mujq) 1 / 2 of this harmonic oscillator. When energies are expressed in units of Tiujq 
and lengths in units of r the normalized trial wave function reads 

(a\ N / 2 SL ( (xt + y 2 k )\ 

where a = a' + za 

In evaluating the expectation value E (Eq. (3)) it is convenient to consider the local 
energy E L (R) = \E'y 1 (_R)i^\E' T (_R) [34]. In general E L (R) is a complex function 

E L (R) = ReE L (R) + UmE L (R) (8) 

with 



N 



ImE L (R) = - a J2((Mn) + eA ext (f k )/c) ■ r k ). (9) 



k=i 



However, evaluation of the expectation value E = J ^t(R) E l (R)^!t{R) dR immediately 
yields 

J ^ T (R) ImE L (R)^ T (R) dR = 0, (10) 

and, therefore, the only quantity to consider in the following is ReE^R). Before proceeding, 
we would like to emphasize that the absolute ground state of the anyon system is a non- 
analytic function of v. Our calculations will simply provide a smooth interpolation. 



III. NON-INTERACTING CASE AND H = 

In the non-interacting case, in the absence of an external magnetic field, the local energy 

is 

-/V 2 2 2 

ReE L (R) = J> + ^±^(1 - a 2 ) + y(A,(f,)) 2 ]. (11) 

k=l 

The expectation value of ReE^R) can be easily calculated for the first two terms of Eq. (11). 

The last term contributes with integrals of the form / ty T (R) —t^tt^^^t(R) dR, which 

fall into one class of N(N — 1) integrals with j = I and a second class of iV(iV — 1)(N — 2) 
integrals with k ^ j,k ^ I and j ^ I. The first class of integrals can be evaluated using [35] 

[°° Ei(ax)e-» x dx = — In - (12) 
Jo n \a J 

with a > 0, Reji > and // > a. (Ei(y) = — j°° y e~ z dz/ z is the exponential integral with 
y < 0.) The result is 
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J* T (R)^* T (R)dRna]n(^ , (13) 



which displays a logarithmic divergence with the cut-off parameter 5 tending to zero. The 
integrals of the second class yield 



/ *t(R) 



Tkj '™ m T {R)dR=-aG , (14) 



\rkj\ 2 \r k i\ 2 



where G = 3 1 / 2 ln(4/3). Putting together all these different contributions one obtains 

E = ^-(m a+-) (15) 



2 V a 
with 

M = 1 + u 2 (N - l)[ln (1) - G(N - 2)] , (16) 
dE 

which attains a minimum (— — = 0) for 

da 

a,=M- 1 ' 2 . (17) 

Thus, the resulting expression for the ground state energy is 

E = N Af 1/2 . (18) 

The logarithmic divergence displayed in Eq when 8 — > has also been found in other 
approximate perturbative treatments of the problem and is widely discussed in the litera- 
ture. To remedy this problem, various solutions were introduced: In Ref. [24] a hard-core 
centrifugal term and in Ref. [25] a pair correlation term were introduced in the trial wave 
function, while in Ref. [26] both modifications of ipr(f) were used. An artificial repulsive 
delta-like potential was assumed in Refs. [27,14,28] when the unperturbated ground state 
wave function is a product of single particle (gaussian) wave functions. Here we assume as 
in Ref. [4] that the cut-off parameter 5 cannot be zero for v > 0, away from the bosonic 
limit, since it corresponds to the square of the nearest distance between the particles. Thus, 
for anyons in the parabolic confining potential 5 is definitely smaller than 1 (in units of Tq). 
In the following we determine 5 by fitting to appropriate results for special values of the 
parameters of the system. 
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FIG. 1. Relative deviation (in percent) of the approximate ground state energy Eq, Eq. (21), 
from the exact ground state energy, E exact , for up to N=72 noninteracting fermions {y = 1) in a 
parabolic confining potential. The dashed-dotted line indicates the asymptotic (N — > oo) value. 



Wu [11] has computed the ground state energy of N anyons in a 2D parabolic potential 
near the bosonic limit v ~ and obtained 



E » [N + N(N -l)u/2). 



(19) 



To regularize the expression for E we make use of this result by expanding E , Eq. (18), 
for v — > and identify the leading term in v 2 with the term linear in v of Eq. (19), with the 
result 



cxp 



1 + vG[N - 2) 



(20) 



With this value of the cut-off parameter the final analytic expression for the ground state 
energy is (see also Ref. [36]) 



E = N[l + u(N - 1)] 



1/2 



(21) 



By construction, it is evident that this formula reproduces the result of Wu [11] in 
the bosonic limit v — > 0. Less trivial, however, is the asymptotics in the fermionic end: 
For large iV it is consistent (up to a numerical factor) with the approximate expression 
E ~ z/^iV 3 / 2 of Chitra and Sen [14] calculated perturbatively from the bosonic end for 
v > 1/N. These authors have also studied the fermionic end v ~ 1 and found for iV ^> 1 
the expression E ~ (8A^ 3 ) 1 / 2 /3. This formula is asymptotic to the exact ground state 
energy E = N c i(l + 8iV c ;) 1//2 /3 for N c i fermions filling the first K closed shells [5], where 
N c i = K(K + l)/2. Note that Eq. (21) provides a monotonically increasing function of v 
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while in the closed-shell case the exact fermionic end has lower energy (by a factor &l 2 fZ) 
than the one calculated from the bosonic end. 

In Fig. 1 we compare exact ground state energies (the sum of occupied harmonic 
oscillator states), for up to iV = 72 fermions (is = 1), with the results obtained from 
Eq. (21). As it turns out, the relative deviation does not exceed 6%. Fig. 2 shows the 
relative deviation for 2,3 and 4 anyons as a function of v (the exact ground state energies 
here are taken from Refs. [1,2,12,13]). In this figure we have considered all the cases for 
which the exact ground state energies are known. When the number of anyons N increases, 
the absolute ground state of the system is a non-analytic function of v because there is 
approximately N 1 ^ 2 number of level crossings [14]. Since our formula for Eq, Eq. (21), has 
the same (N — > oo) asymptotics as the formula obtained by Chitra and Sen [14] one expects 
that, for < is < 1, this relative deviation will be bounded as the number of particles is 
increased. 

It should be noted, that due to our regularization procedure the ground state energy 
obtained, Eq. (21), is not an upper bound to the ground state as one would expect from 
a variational principle. In fact, for N = 2 and v = 1 Eq. (21) yields a value below the 
exact ground state (see Fig. 1). This is a consequence of the fitting of S to the result of Wu 
Ref. [11], which for N = 2 leads to a square root dependence in is, while the exact result for 
this case gives a linear dependence. On the other hand, Eq. (21) applies for the whole range 
of parameters of the system N, is, and u . 




IV. NON-INTERACTING CASE AND H ^ 

In this Section we include an external homogeneous magnetic field. In the presence of 
an external magnetic field H the statistical factor is may change sign because is = e<fi/2irh 
is a fraction of the flux quantum carried by each anyon, 0o — 2nhc/\e\, and this flux 



S 



can be antiparallel to the magnetic field [4,19]. The Hamiltonian is invariant under the 
transformation (xk,Vk, v,(3) — > (xk, —yk, — v, —(3), where (3 = eH/\eH\, and thus the energy 
spectrum is invariant under (u, f3) — > (— u, —j3) (see Ref. [19]). The spectrum only depends 
on \v\, v(3 and the cyclotron frequency to c (apart from N and cu ). 
The real part of the local energy is given by 



N 



ReE L (R) = 

k=i 



a + 



x \ + v\ 



v(3u c 



2uo { 



N 

E 



1 -a 2 + 



0J„ 



4< 



W 



(22) 



We need to compute the contribution coming from the last term in ReE^. To this end, 
we have to solve the integral / ^t(R) ^o ^t(R) dR (N(N — 1) integrals of this form 



contribute to the energy). With the help of the integral [37] 

f°° 2 2 

/ e' a x I y {bx) dx = 
Jo 



IT 1 / 2 A 2 



eSa 2 lv - 

2a 2 \8a 2 , 



(23) 



where Re v > — 1, Re a 2 > (I v (x) is the modified Bessel function), and [38] 



[•CO 

/ e~ ax sinh(fcr) dx = — 
Jo a 



b 2 



with a > b > 0, one gets 



(24) 



N 



*t(R) E rj ^^T(R) dR = 

k,j^k \ rk j\ 

The averaged real part of the local energy is 



-N(N- 1) 



n JV/ f If oo 2 



4< 



u/3u c N{N-l) 



and takes its minimum value for 



do = 1 + 



CO 



2 \ 1/2 



4< 



AT- 1 / 2 . 



(25) 



(26) 



(27) 



The resulting energy minimum is given by (for now we return to standard units of energy 
and length) 



2\ V2 



E \4 • f ) ' Af^-^N(N-l). 



(28) 



As in Section III this expression diverges logarithmically in the limit of a vanishing cut-off 
parameter 5. Following the line of arguments of the previous section, the cut-off parameter 
- representing the squared minimum particle distance - should not be zero except for the 
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bosonic limit v — 0. Having this in mind we determine S by fitting to known exact results 
for the ground state energy. 

To establish these results we calculate the fermion ground state energy from the single- 
particle spectrum of the 2D harmonic oscillator in an external magnetic field perpendicular 
to the 2D plane (Fock-Darwin spectrum [39] ) (see also [33]) 



In Eq. (29) n and m are the radial and angular momentum quantum numbers, respectively. 
The ground state energy of N spinless fermionic particles is the sum of the N lowest single- 
particle energies (Pauli exclusion principle). Following [5] we introduce a parameter z = 
R/P, where R = huo c /2 and P = %(ool + ujI/A) 1 / 2 and express the ground state energy 
in units of P. The parameter z changes between and 1 when the external magnetic 
field is changed between and infinity. The Fock-Darwin spectrum is characterized by level 
crossings. These crossings, occurring at z — z b , z 2 , z^, ...zi, have to be considered in evaluating 
the ground state energy, their number therefore depending upon N. Every interval between 
level crossings is characterized by its own expression for the ground state energy. However, 
only for the intervals < z < z\> and z\ < z < 1 one can write down the expressions for the 
energy as a function of the number of particles N. For the sake of clarity let us consider 
some special values of N. The cases with one and two fermions are not affected by crossings. 
For N = 3 we have one crossing at z b — z x — 1/3 and two expressions for the ground state 
energy: E/P = 5 and E/P = 6 — 3z. This crossing point coincides with the one considered 
in Refs. [5,21] for three anyons. The case with N = 5 has two crossings z b — 1/3 and z\ = 3/5 
and three expressions for energy: E/P — 11 — 2z, E/P = 12 — 5z and E/P = 15 — lOz in 
the intervals < z < Zb, Zb < z < zi and zi < z < 1, respectively. On the basis of these 
special cases one can make the following generalizations: 

• There are N c i — K crossing points z = z b , z 2 , z%, ...z\ for N c i fermions in K closed shells. 
Therefore, there are N c i — K + 1 expressions for the ground state energy of fermions 
and ground state spectra of anyons for this number of particles. 

• One can write the expression Zb = \/(2K — 1) for K closed shells. The last crossing 
point does not depend on K, it is a function of N: z\ — (N — 2)/N and, thus, it is 
applicable for any number N . 

• One can write the expression for the ground state energy in the interval < z < Zb (we 



choose the smallest z b for all particles filling the given shell) or < cu c < u /(K(K — 
1))V2 for K > 2. 



This expression is approximate (to within 6% accuracy (see Section III)) for iV ^ N d , 



where S = (~ N * + 2 -?)> N * is the integer part of [-1 + (1 + 8A^) 1 / 2 ]/2 and 




(29) 



E/P w N 3/2 + zS . 



(30) 



N-Ni-l 



j=0 

N[ = N S (N S + l)/2, and becomes 
closed shells, i.e., iV = N d . 
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• The expression for the ground state energy in the interval z\ < z < 1 or ou c > ujq(N — 
2)/(iV — l) 1 / 2 , determined by the lowest levels from each shell in agreement with the 
Fock-Darwin formula Eq. (29), is 

E/P = N[(N + l)-z(N-l)]/2. (31) 



Having discussed the case of fermions = 1) we now determine the approximate 
expressions of the ground state energy of N non-interacting anyons for the two interesting 
ranges of weak < cu c < uo/(K(K — l)) 1 / 2 and strong cu c > cu (N — 2) / (N — l) 1 / 2 magnetic 
fields. In the weak magnetic field regime it is 



E = PN d [l + \u\(N d 



11/2 



for closed shells and 



E Q = PN[1 + \u\{N - 1)] 1/2 + u/3RS 



(32) 



(33) 



otherwise. Note that these expressions coincide with Eq. (21) in the absence of an external 
magnetic field, i.e., oj c = 0. We find the cut-off parameters by equating these expressions to 
Eq. (28). The result for closed shells is 



5 = - exp 



"A (l + MK^ci " + z 2 /A) - G] + ~rrW + \"\(N d ~ 1)] 1/2 



(34) 



while for open shells 



5 = - exp 



2u(3zT 



u\z 2 T 2 



~ (1 + \v\G{N - 2) + ^^[1 + M (N - 1)] 1/2 + J^L 



(JV-1), 



5 



(35) 



where N > 2 and T = (N - l)/2 + S/JV. 

In the high magnetic field regime (or for weak confinement) the Fock-Darwin single- 
particle energies tend toward Landau levels, the lowest energy state having the quantum 
numbers n = and m < and energy R. In this limit the ground state energy of N 
particles is NR, independent of \v\ [13,20,5]. This exact result can be reproduced with 



M 1/2 = 1 + v(3{N - l)/2 
which gives the cut-off parameter 



exp 



4i> 



[Av{3/\v\ + \v\(N(4G + 1) - (8G + 1))] 



o- 



(36) 



(37) 



Using this choice in the general formula Eq. (28) we arrive at the closed analytic expression 
for the approximate ground state energy 



E = PN[1 + vf3{N - l)/2] - vf3RN(N - l)/2. 



(38) 



Besides the high magnetic field (or weak confinement) limit used here to fix the cut- 
off parameter 5, the expression obtained for E reproduces the exact ground state energy 



11 



of non-interacting fermions (the sum of the TV lowest Fock-Darwin energies) in the whole 
magnetic field range beyond the last crossing of the iVth level with n = 0, m = — (N — 1), 
which defines the so-called maximum density droplet [7]. Finally, approximate expressions 
Eqs. (32), (33) and (38) give the ground state energy E = PN of iV bosons {y = 0) in a 
magnetic field and harmonic confining potential. 

The separate discussion provided here for small and large magnetic fields is in correspon- 
dence with the treatment of the N = 3 case discussed in Refs. [21,5]. 



V. COULOMB-INTERACTING CASE AND H = 



We now include the effect of the Coulomb repulsions between anyons 



ro 



N 

E 



m 



2ob fcJ#fc |r fc j| 

the expression for the real part of local energy ReE L (R), Eq. (11), but in a vanishing 
external magnetic field. Here, as in Section III, we assume that v = \v\. 

The Coulomb interaction part contributes with N(N — 1) integrals of the form 

J ^t(R) t^—;^t(R) dR. These integrals can be evaluated using Eq. (23) and [35] 



\r kj 



f 

Jo 



~I u {bx) dx 



b u 



Va 2 - b 2 (a + ^a 2 -b 2 Y' 
where Re v > — 1 and Re a > \Re b\. The result is 



(39) 



2 J 



[ ^ t {R)t^-^t{R) dR - 

J \ r ij\ 

The averaged (real part of the) local energy is 

E=^(U a+- + 2M a 1 ' 2 
2 V a 



nay/ 2 



(40) 



with 



M 



7r\V 2 N- 1 r 
2/ 2 ~a~ B 



(41) 



(42) 



dE 



The extremum condition — — = leads to the equation 

da 

X 4 - MX - N = 



(43) 



for X = I /a 1 / 2 . Two complex and two real solutions of this equation can be found by the 
Descartes-Euler method [40]. The minimum energy is given by the expression 



En — 



N 



'2 




2M 



(44) 



and it is achieved at the point 



12 



X = (A + B) 1 ' 2 + [-(A + B) + 2(A 2 — AB + B 2 ) 1 ' 2 } 1 ' 2 



(45) 



where 



A 
B 



M 2 /128 + ((AT/12) 3 + (M 2 /128) 2 ) 1/2 
M 2 /128 - ((AT/12) 3 + (M 2 /I28) 2 f /2 



1/3 
1/3 



(46) 



Again, the ground state energy E , Eq. (44), has a logarithmic divergence in the limit 
8 — > 0. Assuming that A/" can be regularized, one can recognize two limits of interest. One 
corresponding to weak correlations, r /a_B "C 1, 



En^N M l/2 + 



A/" 1 / 4 

and another for strong correlations, ro/as 1, 



(47) 



(48) 



In order to determine the cut-off parameter 5, and due to the lack of analytic results, we 
need to fit to known numerical results for the ground state energy at special values of the 
parameter tq/qb- 
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r /a B 





N = 10I 




A/ = 9| 




AT = 8| . 




N = 7\ . 



FIG. 3. Coulomb interaction parameter ro/as dependence of the ground state energy for 7 - 
10 electrons calculated by variational [41] and fixed-node quantum Monte Carlo methods [42] - 
dashed curves (results of both calculations are indistinguishable in these curves) and by formula 
(44) - solid curves. 
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In Fig. 3 we compare the ground state energies calculated for 7-10 electrons using Eq. 
(44), with the non-interacting M = l + v(N — 1), to variational [41] and fixed-node quantum 
Monte Carlo calculations [42] (see also [43]). 

From Eqs. (47) and (48) follows that the contribution of the statistics (dependence upon 
v) in the ground state energy is important for weak, and negligible for strong, Coulomb 
correlations to/ob. For large values of r^/as one can compare the dependence of the ground 
state energy, Eq. (48), with the estimate given in Ref. [44]. The asymptotic behavior of the 
ground state energy with N derived from our expression is E ~ iV 5 / 3 (as in Ref. [44]). 



VI. COULOMB-INTERACTING CASE AND H ^ 



Finally, we consider the case of a confined Coulomb anyon gas in an external magnetic 
field. The resulting expression for the averaged local energy is 



E 



N 



Af 



a 



1 L u c 
- i + 



a 



4^ 



: \2M a «A-^ N[N - l) 



2uj 



(49) 



Following steps similar to previous sections we obtain the ground state energy (in standard 
units) 



En 



N 



2u>o Xq 



(50) 



where Xo is formally the same expression as Eq. (45) after replacing M — > M\ = A/"/(l + 
lv^/(4uj^)) and M — > Mi = M / (I + / (4uJq)) in Eq. (46). The asymptotic expressions are 



Eq NhuJc 



Ml 12 



and 



Eq 



3NhuJ t 



o 



2/3 



1+ ^ 



1 + 



u(3uj c (N-l)' 
2uj 2 



vPu c (N-l) 
3u 2 



(51) 



(52) 



for very small and very large values of ro/a^, respectively. We regularize the logarithmic 
divergence by fitting 5 to known numerical results for the ground state energy of quantum 
dots in external magnetic fields. 
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FIG. 4. Magnetic field, H, dependence of the ground state energy for N = 3 and N = 4 
spin-polarized electrons in a harmonic potential calculated in Ref. [45] (the dashed curves), and 
using Eq. (50) (the solid curves). As in Ref. [45] we used fkvo = 3.37 meV (ro/ae = \JH*/ {Kuq), 
where the effective Hartree H* is equal to H* ~ 11.86 meV). 

We compare our results with the ground state calculations of Ref. [45] for GaAs dots 
in an external magnetic field. The results of these calculations are very close to the results 
of Ref. [46] and Ref. [47] computed by exact diagonalization and quantum Monte Carlo 
methods, respectively. For GaAs quantum dots M* = 0.067M, and the dielectric constant 
is e = 12.4. Therefore, the effective Bohr radius a* B = h 2 e/M*e 2 is a* B ~ 97.90Aand the unit 
of energy - the effective Hartree (H* = M*e 4 /(e 2 ft 2 )) is H* ~ 11.86 meV. The cyclotron 
frequency is u c = eH/(M*c) (for simplicity, in this part of our work we assume v = \u\, 
e = \e\ and H = \H\ while a correct combination of the signs of these quantities is given 
in Eqs. (50) - (52)). Thus, the energy quanta for this frequency is hu c = 1.7269 • H- 
(meV/T). Here we took into account that Tiu c = heH/(M*c) = 2M/j, B H/M*, the Bohr 
magneton is fiB = eh/(2Mc) = 0.05785 meV/T and the magnetic field H is measured in 
Tesla (T) magnetic units. The Coulomb interaction parameter r^jas in our case is equal to 
r /a B = ^H*/(hLj Q ). 

To compare with the results of Ref. [45] for spin polarized electrons, we calculated the 
ground state energy using Eq. (50) with the expression M — N (i.e., non-interacting M 
with v = 1), for three and four particles as a function of the magnetic field strength H . As 
in Ref. [45] we considered %Uq = 3.37 meV. Comparison of these results is displayed in Fig. 
4. The deviation of our results with respect to the ones given in Ref. [45] is no more than 
10%. 

It turns out that this expression for M with v = 1 is appropriate for the description of 
a small number of electrons (N = 3, 4, 5, 6) and not large magnetic fields. For large number 
of particles and a wide range of magnetic fields one can write the approximate expression 
for M (here we return to the original signs of u, e and H) 
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M = F 
and thus for 5 

1 



5 = - exp 



fl + M G(iV-2) + M^ 
V 4w 2 



with 



'vPuJcN 1 ' 2 ' 



1/21 



W 



(53) 



vf3uj c C 
2\v\ujq 



Q 



(54) 



1 + 



2, ,2 ' 



(55) 



iV 2 - 2iV - 3 
4(iV- 1) 



+ \u\ 2 G(N-2), 



(56) 



2ATV2 

c = ]7ri + (i + M(^-i)) 



1/2 



UJQ 



1/2 



(57) 



and 



Q = 2 



v{3uj c N 1 I 2 

\v\ 2 0Jq{N — l) 5 



;i + lH(^-i)) 



1/2 



(58) 



for > 2. Here we took into account that A/" depends weakly on the Coulomb parameter 
r /a B (the results indicated in Fig. 3 have been calculated with M not having this parameter 
dependence) . 

In Fig. 5 we compare the ground state energy calculated with the expression Eq. (50), 
(using Eq. (53) for M with \u\ = 1) for 16 < iV < 40 and < cu c /cj < 20 (r /a B = 1.911), 
with the calculations for a classical system of electrons of Ref. [48]. The deviation is maximal 
(no more than 15%) in the range l/(K(K - l)) 1 / 2 < u; c /u; < (N - 2)/{N - l) 1 / 2 , where K 
is the number of closed shells. This range of magnetic fields corresponds to the crossings of 
Fock-Darwin levels (see Ref. [39] and Section IV) and, therefore, the single particle ground 
state energy of electrons changes many times as uo c increases. One can suppose that the 
Coulomb interaction shifts the levels but the qualitative structure of the many particle 
ground state is still complex. Thus, in this range of parameter uj c / ojq the expression for M 
is not uniquely defined. We could not find a more appropriate expression for J\f than Eq. 
(53), for the magnetic fields indicated in Fig. 5. 
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FIG. 5. Ground state energy = (Eq — Nfko) / (Hujq) for 16 - 40 electrons calculated using 
the expression Eq. (50) for ro/a# = 1.911, applying the expression for M Eq. (53) with \v\ = 1 - 
solid curves, and energy for classical electrons (Ref. [48]) - dashed lines. Here u = (u>q +lv1/A) 1 / 2 
and t = u> c /u>o. 



VII. CONCLUSION 

We have used the anyon concept combined with a variational calculation to obtain an 
analytic closed-form expression for the approximate ground state energy of N non-interacting 
and Coulomb-interacting particles in a 2D harmonic confining potential, with and without an 
external magnetic field. The crucial point of this approach is the appearance of a logarithmic 
divergence connected with a cut-off parameter, when evaluating the contribution of the gauge 
field vector potential. Following arguments from the literature, according to which the cut- 
off parameter cannot be zero (except for the bosonic limit), we used it to fit our results to 
exact and numerical ground state energies known for special values of the system parameters. 
In doing so we provided closed analytic expressions for the approximate ground state energy 
depending upon N,uo,r /aB and uj c /ujq. 
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